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Abstract. We consider cosmological solution for Einstein gravity with massive fermions with 
a four-fermion coupling, which emerges from the Hoist action and is related to the Barbero- 
lmmirzi (BI) parameter. This gravitational action is an important object of investigation 
in a non-perturbative formalism of quantum gravity. We study the equation of motion for 
the Dirac field within the standard Friedman-Robertson- Walker (FRW) metric. Finally, we 
show the theory with BI parameter and minimally coupling Dirac field, in the zero mass 
limit, is equivalent to an additional term which looks like a perfect fluid with the equation of 
state p = wp, with w = 1 which is independent of the BI parameter. The existence of mass 
imposes a variable w, which creates either an inflationary phase with w = —1, or assumes an 
ultra hard equation of states w = 1 for very early universe. Both phases relax to a pressure 
less fluid w = for late universe (corresponding to the limit m — > oo). 
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1 Introduction 

The investigation of Dirac fields in curved spacetime has been addressed in several works for 
many years, especially in the last decade (see, for example, Refs. [1-10]). There are many 
papers considering not only the classical aspects of the theory, but also quantum ones. In the 
present work, we deal with the classical aspects of massless Dirac particles in the spatially flat 
Friedmann-Robertson- Walker background in the presence of the Barbero-Immirzi parameter 
[11, 12], refers to as BI parameter in the following. This parameter, described by the Hoist 
action term [13], was introduced under the non-perturbative quantum gravity perspective, 
and it represents then a new dimensionless parameter coming from a more fundamental 
theory. 

The vacuum gravitational action together with Hoist term can be written as 



where k = 16ttG, e a p^ u is the Levi-Civita tensor and /3 is the BI parameter, which should be 
real and positive. In General Relativity (GR) this parameter makes no effect on dynamical 
equations, because e a p flu R a P f11 ' vanishes due to the cyclic symmetry of the Riemann tensor, 
-gafijiv _|_ ftavPfj, _|_ fta/iup = Q However, this is not so if torsion is present. The simplest way to 
introduce torsion is to consider the Einstein-Cartan action coupled to a Dirac field, such that 
the equations of motion for torsion become non-trivial. For an introduction and review on the 
theories with torsion, one can see [14-16], the last reference gives nice overview of the relevant 
cosmological models. It is worth noticing that, in this scenario, where fermionic matter is 
present the BI parameter and/or torsion can affect the gravitational dynamics, providing an 
interesting way to investigate its classical and quantum effects. It was first proposed in Ref. 
[17] that the Immirzi parameter leads to an effective four-fermion interaction mediated by 
the Immirzi parameter. The same effect can be achieved in the framework of Einstein-Cartan 
theory with minimal or nonminimal coupling of fermions to torsion. 

The first studies of the effect of the Immirzi parameter was done using the minimal 
coupling procedure (MCP) of the fermions to gravity [17, 18]. However it was pointed out 
first in [18] that the effective four-fermions interaction depends on the choice of the coupling 
of fermions to gravity and that we can assume non minimal coupling procedure. 
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This is due to the fact that the action describing Dirac theory in flat spacetime is 
invariant by the introduction of a total divergence, say, d^V^ 1 with = ip^^ip or = 
ip'j^jsip. However when torsion is present, there may be an extra coupling, for such terms 
produces extra term T^V^ in the action (here = is the trace of the torsion tensor). It 
means that a whole class of equivalent actions in flat spacetime are related (though MCP) to 
inequivalent actions in curved spacetime. This problem was especially addressed by Kibble 
[19] (see also [15]). 

In [18] a one complex parameter family of non minimal coupling was investigated, in 
[20] another family of non-minimal coupling was proposed while in [21] a general 3 parameter 
family of non minimal coupling was investigated. In our note we assume that the fermion 
coupling is the minimal coupling and we investigate the consequences of this coupling on the 
cosmological evolution. 

It can be argued that the minimal coupling is the most natural coupling of gravity to 
fermion based on the fact that non minimal couplings are sourced by components of the 
torsion that do not appear naturally in models of spinning matter. For instance it is possible 
to treat fermions, at a large scale, as a fluid with an intrinsic spin density. Following the 
Lagrangian formalism for this exotic spin fluid [22-24], one can see that the contribution 
from such fluid always produces a traceless torsion [25] (for the extensive discussion of the 
theories with spin fluid in cosmology, see also Refs. [26-32]). In this case the trace of the 
torsion tensor = T°^ a has smaller physical meaning and its presence is completely due to 
the nonminimal procedure [14, 15]. It turns out that the non-minimal procedure [33] 



(we use the notations of [15], correspondingly is the Riemannian covariant derivative, 
without torsion) is the most consistent one, especially at the quantum level. Of course it is 
possible to view this procedure as a MCP with a modified connection [34]. In any case, the 
result will be the equivalent action with current-current interaction found in [17]. For the 
corresponding discussion of the issue related to the BI parameter effects, see, for example, 



In the present paper we consider a particular cosmological solution, with a FRW-like 
metric, within the standard approach to BI parameter, and show that this parameter does 
affect the cosmological solution. 

The paper is organized as follows. In the next section we briefly discuss the existing 
ambiguity of the MCP in Einstein-Cartan theory. In Sect. 3 the details of the derivation of 
equations for the metric and for the components of the Dirac field are given. The details of 
the derivation of the Energy- Momentum Tensor are settled in the Appendix A. In Sect. 4 
we discuss the cosmological solution and verify the consistency conditions. The calculations 
are not completely trivial and hence we present them in some details. Finally, in Sect. 5 
we draw our conclusions and discuss possible perspectives of cosmology based on the gravity 
theory with BI parameter. 

2 Derivation of dynamical equations 

The action for gravity theory with an additional Hoist term with the BI parameter and 
minimal fermion coupling, in the presence of massive Dirac fields can be written in terms of 




(1.2) 



Refs. [34-36]. 
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a spin-spin coupling in the form [17, 18] 



S = Jd^xy/^f—R+^^V^-V^i/i) -m^-ej^X , (2.1) 



where 9 is the spin-spin parameter related to Immirzi parameter j3 by 

32 /3 2 + 1 

and the axial current is defined as 



(2.2) 



= ^r 57 ^v • (2.3) 

Here, T§ = iroril^ra is the chiral Dirac matrix. The Dirac matrices, 7^, are expressed 
in terms of the vierbeins e a v and the Dirac matrices T a in tangent flat space according to 

Two important observations are in order. The same action (2.1) can be seen as the 
result of integrating out completely antisymmetric torsion in the Einstein-Cartan gravity 
coupled to the fermion fields forming axial current. In this sense all our further consideration 
can be attributed to one of the two models, namely Hoist gravity with BI parameter and 
Einstein-Cartan gravity. Further aspects of the difference between these two formulations 
were discussed in [20]. 

The second observation concerns the relevance of the choice of the MCP procedure, 
which we have already mentioned in the Introduction. The MCP corresponds to the value 
771 = 1/8 in (1.2). Taking into account the relation of our model with the Einstein-Cartan 
gravity, it becomes obvious that the transition to an arbitrary non-minimal interaction be- 
tween fermion and torsion is performed by the multiplication of 9 in (2.2) by a factor of 64r]f. 
Therefore the "nonminimal" value of the parameter should be 

^nonminimal — ^ _|_ j ' (2-4) 

In the ideal world with only one kind of non- interacting fermions, this would be the end 
of the story, because the single parameter 9 contains the information about both /3 and 
r/i. However, in many physically interesting theories such as Standard Model of particle 
physics, there are different kinds of fermions, all of them interacting to scalar (Higgs) by 
means of Yukawa interaction. At quantum level the nonminimal parameters 771 for different 
fermions are renormalized, become running parameters, and their running depends on the 
values of Yukawa couplings [33]. As a consequence, the values of rj\ must be different for 
different fermions. In what follows we present a simplified consideration with only one 9, 
corresponding to the one single type of fermionic axial current J^. One can regard this 
approach as simplification, which is based on the assumption that the running of different 
parameters r/i is not very strong and that MCP hypothesis can serve, after all, as a good 
approximation. 

The variation with respect to the Dirac fields gives the dynamical equations 

vfV^ = WJ^T^if) + , (2.5) 
iV M $y = - 29J^T 5 ^ - w4 ■ (2.6) 
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The Einstein equations, coming from the variation with respect to g a s, can be written 

as 1 

1 k i - 1 

R a p - 2 R 9aP = G afi = = ^ « [^7(«^)^ - V^T^] - -Kg af3 C , (2.7) 

where 

C= l - (ff V M V - V^fty) - m^V - , (2.8) 

and 7( Q Vg) means \(^f a V p +7,gV a ). By using equation (2.5), we obtain 

i -9 
G a /3 = [^T(qV^)^ - V (/ #7 q) V>] - -KgapJ^JV . (2.9) 

Let us note that our expression for the energy-momentum tensor (2.7) is slightly different 
from the one of [1], which was obtained by variation with respect to the vierbein. The 
expression (2.7) has been obtained by variating with respect to the metric. The technical 
details can be found in the Appendix A. 

The next step is to calculate the equations for the specific metric corresponding to the 
homogeneous and isotropic universe (FRW space-time). The spin connection, 

<A = \ (e ba d^e a a - e aa d^e b a + e aa d a e\ - e ba d a e\+ 

+ e bu e aX e cll d x e c u - e av e bX e Cfl d x e c v ) , (2.10) 

can be obtained for the FRW spacetime with null spatial curvature, ds 2 = dt 2 — a 2 (t)(dx 2 + 
dy 2 + dz 2 ) = a 2 (r])(dr] 2 — dx 2 — dy 2 — dz 2 ), by taking into account the vierbeins 

e° = (afa), 0,0,0), e a ! = (0, a(r?), 0, 0) , e\ = (0, 0, a(r?), 0) , e a 3 = (0, 0, 0, a( V )p.ll) 

Here we consider the conformal time, rj, defined by a(rj)drj = dt. Thus, the non-null spin 
connection components ui ab fl are 

,,01 _,,02 _ ,03 _ a ' ,,10 _ , ,20 _ , ,30 _ a ' ( -, 9 \ 

UJ i — UJ 2—U 3 — 7T~ > W 1— W 2 — W 3 — -7— , {2.12) 

2a 2a 

where a' = da/drj. 

3 Cosmological solution and consistency conditions 
3.1 Dirac equation 

In order to solve the equation in a cosmological background, we make the assumption that 
the fermions are spatially constant, diip = 0. With this assumption and using the form of 
the connection (2.12) we get 

a' a' - 

V ^ = 9 ^, Vo^ = 3 V, = — IW , = —4>T Ti , i = 1,2,3. (3.1) 

2a 2a 



1 The expression for the Energy-Momentum Tensor for the free case is known from [37]. 
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The Dirac equations (2.5) and (2.6) can be written as 



-r (V + ^V) = m^ + 26{Y 5 T a ^)J a 
a \ za / 

--(& + -^-4>Vo = m^ + 26 J a (i>T 5 T a ) . (3.2) 
a \ la J 



where J a = (ipT 5 T a ip). 

From these equation we can extract the conservation law: 

$r V)' = -— (W). (3.3) 

a 

which is identical to the current conservation 

o = v M $YV) = d^ri?) + rv^yty) > (3-4) 

where T^p^ is the Riemannian affine connection. This can be easily integrated out 

M 3 

V (V) = — 
cr 

where M is a mass scale that gives the size of the condensate and V a = ijjTaip is the vectorial 
current. 

One can also get the following equations for the symplectic current: 

^ (VToV' - ^ToVO = 29 J - r + m W> (3-5) 
^ (#^' - $T^) = 0. (3.6) 

In order to analyze further the spinor equations of motion, let us introduce not only the 
vectorial and axial current V a and J a but also the pseudo scalar P = i^T^ip, scalar S = ipip 
and Lorentz tensor L a & = ^\Tai ^b]tp currents. 

We can check that the equation of motion for the scalar 5 and pseudo scalar P and 
the time component of the axial current form a closed subset of equation. Indeed, in the 
Appendix B we derived the closed set of equations 

^(a 3 S)' = 2a 3 9PJ , (3.7) 
la 

^(a 3 J )' = -a 3 mP, (3.8) 

— (a 3 P)' = a 3 {m J - 29SJ ) . (3.9) 

Zi(X 

From these equation we can conclude that a 6 (S 2 + J% + P 2 ) is a conserved quantity. Moreover 
we have the continuity equation: 

[a 3 (0J o 2 + mS)]' + 3a'a 2 9J$ = 0. (3.10) 

We can also show that the spatial components of the axial vector are conserved, 

(a 3 J 4 )' = 0. (3.11) 
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For completeness we can also give the evolution equation for the other currents: 



(a 3 Vi)' = mL 0i + 29e iaj J a V j , 

1 

± (a 3 !/)' = 29{e^J,L k - 2J L 0i ) , (3.12) 



(a 3 L 0i )' = -mVi + 26e iab J a L oh - 6J e i]k V k , 



where L l = l/2e^ k Lj k , and e a ij is the 3-dimensional Levi-Civita tensor (e a ij = if a = 0). 
One can extract from these equations one conserved quantity, namely 



a^W + L^Loi + UV). (3.13) 



3.2 Einstein equation 



The spatial and time-space off-diagonal components of the fermion energy momentum tensor 
are given by 

T H = \ (^7(iV i} V - V (i ^ 7i) V) + a 2 % -0J% = , ijtj, (3.14) 

1 — — 

Toi = - (^7(o v i)^ - V (o^7i)V>) , (3.15) 

Too = \ (^ToVqV - V ^7o^) - a 2 0J% . (3.16) 

Direct algebraic manipulation with Eq. (3.14), using the derivative (3.1) together with the 
identity r^r^To — Tor^r^ = lead to the conclusion that the first term in (3.14) vanishes. 
One can also easily see that the condition (3.15) is proportional to the LHS of (3.6) and thus 
also vanish. Finally, thanks to (3.5) and (3.1) one can see that the first term of (3.16) is 
equal to a 2 (29J a J a + mtpi/j). 

In summary this shows that one can express the energy-momentum tensor of the fermion 
in the standard perfect-fluid form, 

Tpv = (p + p)u^,u u - pg^ v , (3.17) 

where, for the FRW metric, = (a, 0, 0, 0). Using identifications Too = po 2 an d Tij = 
paPrjij, we are able to express p and p as 2 

p = 6J 2 + mS and p = 9 J 2 . (3.18) 

where J 2 = J fl J» = J a J a . 

The energy conservation law reads 

p + 3H (p + p) = with H = a/a, (3.19) 

This is consistent with the Dirac equations of motion derived previously. 



2 Let us note that the zero Energy-Momentum tensor without 8 in the massless case is only due to the fact 
we perform this calculation for a conformally-fiat metric. In this case the trace of the Energy-Momentum 
tensor completely controls all its components. Some extended discussion of this issue was recently given in 
[38]. 
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For the FRW spacetime, every off-diagonal component of Einstein tensor vanishes iden- 
tically, Goi = Gij = (i 7^ j). This is consistent with the form of the energy-momentum 
tensor just calculated above. The temporal and spatial components of Einstein equations 
(2.9) read, respectively, 

^a 2 {9J 2 + mS) = \a 2 p, (3.20) 
VflJ 2 = Vp, . (3.21) 

while taking the trace of the Einstein equation, we can express the acceleration in terms of 
the physical time coordinate t as 

Q A = _^ A6 J 2 + mS) . (3.22) 

3.3 Massless case 

In the case where m = we can see that the BI parameter for a massless Dirac field produces 
the same effect as of a perfect fluid characterized by the equation of state (EOS) p = p = 9 J 2 . 
It is remarkable that, if m = 0, both energy density and pressure in eq. (3.18) do depend on 
9, while their ratio does not. Hence the theory under consideration does not have a smooth 
limit for 9 — > 0. In this case the energy conservation law (3.19) gives us a new scaling law 
for energy density of massless fermions in a gravity theory with non-zero BI parameter, 

p~a" 6 . (3.23) 

In the considerations presented above we did not take into account the energy density 
and pressure of the free part of the fermion spinor field. According to our calculations this 
part is zero. However, this apparent result is only due to our specific treatment of massless 
fermionic field in the cosmological setting, where we are essentially looking for the non- 
conformal part of T^ v . The conformal part which does not couple to the conformal factor 
a(t) can not be seen in this approach for the following reason. We were looking for the time- 
dependent spinor field, but this is definitely not a right idea for the massless field which can 
not be done space-independent by any choice of the local reference frame. In order to see this, 
in the conformal case, it is sufficient to consider the flat-space limit. The space-independent 
fermion satisfies the equation 

r°«9o^ = 0, (3.24) 

that simply means doip = and constant field. Then the Eqs. (3.14), (3.15) and (3.16) 
tell us that Tn V = 0. However the right way to take the energy density and pressure in the 
conformal part of fermion T^ v into account is quite different. One has to remember that 
in the massless field case the dependence on time and space coordinates should be related, 
because the velocity is universally fixed for all modes of the field. For example, the free 
solution of iT^d^ip = can be always considered as a superposition of massless plane waves 

ll>(x°, X*) = V e *( fe 0*°-kr) ^ k 2 _ k 2 = q _ (3 25) 

It is an elementary exercise to check that each one of these waves creates pressure which 
is 1/3 of its energy density. Therefore the total pressure po and energy density po of the 



3a 



12 



a a 2 
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free fermionic field are related as po = po/3. Taking into account our previous result we 
arrive at the following energy density and pressure of the fermion field in the presence of BI 
parameter: 

p = p + 9 J 2 and p = Po + 9J 2 = - p + 6 J 2 . (3.26) 

Let us remark that the energy density and pressure of the free radiation content of the Uni- 
verse is not necessary fermionic one, but can also include other fields, such as electromagnetic 
radiation. The unique assumption which was done here is that these two quantities are re- 
lated by the equation of state po = po/3, and this is indeed satisfied for the asymptotically 
free ultra-relativistic fields independent on their spin. 

The Eqs. (3.26) indicate that the matter content of the early Universe in the presence of 
the Hoist term is characterized by a sum of two fluids, one of them is conventional radiation 
and another has equation of state which is independent on the value of rj. Of course, the 
relevance of this second component depends on 9, but the equation of state does not. In 
order to have better understanding of the new terms with 9 J 2 in (3.26), we can consider the 
simplified cosmological model with vanishing po and pq. It is easy to solve the Friedmann 
equation, H 2 = (8irG/3) p and obtain the corresponding rule for the expansion of the 
universe, 

a ~ t 1/3 , (3.27) 

which, again, is different from the laws of expansion for free radiation (with a ~ t 1 / 2 ) and dust 
(with a ~ i 2//3 ). The results (3.23) and (3.27) do not depend on the value of BI parameter. 
The only one important requirement is that this parameter should be nonzero. However, 
the scaling law (3.23) shows that the new terms decay much faster than radiation and hence 
become irrelevant in the course of expansion of the Universe. At the same time, there is a 
chance to see the traces of such term in observational cosmological data, especially in the 
ones related to the cosmic perturbations. In order to elaborate this idea one has to develop 
the complete cosmological model with BI parameter. We postpone this issue for the future 
work. 

3.4 Massive case 

In the massive case one has to distinguish the early and late time regime. Let us recall that 
according to the Dirac equations we have the conservation rule 

M 6 M 6 

where M, M are mass scales that control the size of the fermion condensate. This shows that 
S is necessarily decreasing as the universe expand, The late regime of the system is attained 
when 

o 29M 3 

a 6 >> . 

m 

In this regime we have that 29S <C m which means that we can neglect the non linear term 
in Eq. (3.9) and the equations for the evolution of Jo, S, P are dominated by the mass term; 
the influence of the 9 coupling can be ignored. We have, in this case, 

d t {a 3 S)^0, (3.29) 
d t (a 3 J ) ps -2m a 3 P, (3.30) 
d t {a 3 P) » 2m a 3 , J . (3.31) 



J 2 + P 2 + S 2 = — , J? = — (3.28) 
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This system of equation can be easily solved to give 



M 3 M 3 M 3 

J = — L cos [2 m (t - t )} , P = — ± sin [2m(t - t )] , S = — £ , (3.32) 

or a 13 a J 

where Mi and M2 are mass scales characterizing the condensate and satisfying M 6 = Mf + 
M| . In this late time evolution the condensate evolve like a pressure less fluid with p = 
where the energy density is dominated by the scalar component and evolve like p ~ a~ 3 while 
a-t 2 ' 3 . 

Another way to understand this result comes from the fact that the energy density 
p = 9 J 2 + mS possess two components. The scalar components scales like S ~ a -3 while 
the current component scales like J 2 ~ a -6 therefore at late time the energy contribution is 
dominated by the scalar component. 

On the other hand this means that at an earlier time the current is going to dominate 
the dynamic of the cosmological evolution. The crossover time takes place when we can no 
longer neglect the influence of the non linear 9 term in the fermionic evolution equation (3.9), 
that is when 

, 26M 3 

a 3 . 3.33 

m 

At this crossover time both mS and 9 J 2 are of the same order of magnitude provided we 
assume that M%,M2 and M are of comparable magnitude. The energy density at this time 
is of the order m 2 /9 that is of the order (/3 2 + 1)/ f3 2 m 2 Mp where Mp is the Planck mass. 

There are then two radically different early time evolution depending on whether at 
this cross over time the current is space like or timelike. Let us assume first that at this 
cross-over time the current is timelike: i-e J5 3 — J 2 > and that 9 J 2 is of the same order as 
mS. In this case the energy density is dominated at an earlier time by the current term 9 J 2 , 
the fermionic field becomes effectively massless and the dynamic for the fermionic equation 
of motion is entirely dominated by the non linear term and given by the massless equations: 

^(a 3 S) f w 29J (a 3 P) , 
la 

^(a 3 J )'«0, 

— (a 3 P)' w -29J (a 3 S). (3.34) 
2a 



In this regime applicable for early time 



(3.35) 



m 



the solution reads 



a;' 





where M 6 = Mf + M|. In this early time regime the condensate evolves as massless conden- 
sate p = p with p ~ a~ b while a ~ t$ . It is quite remarkable that the early time evolution 
although dominated by the presence of 9 is following an equation of state independent of it. 
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Finally, there is another early time regime accessible to our system. This regime happens 
if J 2 < 0, i.e., the current is space-like at the cross-over time. In this case we cannot have that 
the current term 9 J 2 dominates the scalar contribution since we should respect the constraint 
that p = 9 J 2 + mS > 0. What happens in this case is that the quantity p + p = 26 J 2 + mS 
which is positive at late time decrease to eventually vanish. When this quantity vanishes, 
and since p = —3H(p + p), the energy density becomes constant. This means that we enter 
an inflationary phase provided the value of the density energy is not zero. This shows that 
the other regime of the theory, characterized by a spacelike current at early time correspond 
to an inflationary era that relax to pressure less dust at late time. The value of the effective 
cosmological constant at early time depends on the initial condition but it can be bounded 
by the value of mS at crossover over time. Indeed, since the inflationary era is characterized 
by 26J 2 = —mS, the effective cosmological constant is p\ = mS/2. The value of mS/2 at 
early time is necessarily smaller than its value at the cross over time which is m 2 /29 2 . So 

2 

PA < ^2- ( 3 - 37 ) 



4 Conclusions 



We have considered the cosmological solution for the metric-spinor gravity with the Hoist 
term, especially the effect of the non-zero BI in a cosmological setting and shown that there 
are FRW-compatible solutions. One of the most remarkable result is that the EOS of the 
self- interacting spinor matter does not depend on the BI parameter in the massless limit, 
if we disregard the effect of free massless fermions. For the massive case we have identified 
two different regime. In the first regime the theory is effectively massless at early time and 
behaves as a perfect fluid with an ultra-hard equation of state w = 1. In the second regime the 
fermionic matter behaves effectively as a cosmological constant and creates an inflationary 
phase which is relaxed at late time into a pressureless fluid. 

It would be definitely interesting to check whether the same effects take place for the 
spinning fluid. In the positive case this may lead to a potentially observable consequences for 
the early Universe. The effect of torsion and/or Hoist term on the EOS for the hot matter 
is depending on the existence of axial current and on an arbitrary parameter defined in 
Eqs. (2.2) or (2.4). In principle, some cosmological observations can be helpful in getting an 
upper bound for 9. Let us note that recently the discussion of the effect of Hoist term and/or 
torsion on the difference between EOS for photons and hot fermions (quarks and leptons) 
has been discussed in [41] in the framework of free fermion theory and loop quantum gravity. 
The motivation for this study was to see whether the fine balance required by Big Bang 
nucleosynthesis (BBN) holds in the presence of loop quantum gravity effects. It was shown 
that the possible violation of such a fine balance due to the possible effects of torsion or 
Hoist term is very small. In our understanding, the real physical situation can be even more 
simple, because fermions and photons are supposed to be interacting with each other at the 
BBN epoch. Such an interaction can control the balance between fermions and radiation 
through emission and absorbtion of photons by fermions and also by quantum effects such as 
creation and annihilation of fermion-antifermion pairs by radiation. As a result, the EOS for 
the matter content of the Universe should be treated as unique at that epoch and the effect 
of loop quantum gravity or torsion can not probably change, in principle, the fine balance 
between expansion rates of fermions and radiation. At the same time, this total EOS can 
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be affected by the mentioned manifestations of a new physics, including the four-fermion 
interaction. 

Note added. When we were preparing this manuscript for submission, the preprint with 
partially similar content [40] has been published. 
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Appendix A. Derivation of the Energy-Momentum Tensor 

For this end one has to take 

9 fi u^9'=9 fi u + h flu , ^ = gr-bf v , V /= ?=v / =5fl + 5), (4.1) 
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where only the contributions up to the first order in hF" were kept, and h = h^a = h^g^ 11 . 
The expansion for the vierbeins are [39] 

e' ba = e ba - \h a p e b ^ , e'\ = e\ + \hje a p , (4.2) 

such that J'^J'^ = J^J a . This means that 5{JPJ p )/5g^ = 0. 

In order to find the energy-momentum tensor coming from the kinetic part of the Dirac 
action 

S k = % - J <fx^C k = l -j d^x^-g {^fV^ " V^W} , (4.3) 
we shall consider its variation in h^w By straightforward calculations we obtain 

~g(} + \){^ + \h p ui P )d^{g^ - hn 

Ma b (lv + \h\i p )^ v - h^)} , (4-4) 



S' k = \ d* Xy 



-1- , ,' ab 

where we used the covariant derivative of Dirac spinors in terms of the spin connection, uj ab fl , 

V M V = + l -u a \a ab ^ , = 8^ - l -to ab ^a ab , (4.5) 

and (T a b = \ [r a , Tf,} . Collecting the terms up to first order in h pu , and expressing the spin 
connection as 

, ,'ab _ ab , a ab 



- 11 - 



we find 



+ 1 / d 4 x^ {fy"<^ + VW^} A^ ab M . (4.6) 



It is easy to verify that Aw"^ can be written as 



Aw«*, = i {v^ (e aa e b ^h a ^ - (e ba e^h a ^ } . (4.7) 

Replacing this expression into equation (4.6) and performing integration by parts, one can 
conclude that after some algebra all contributions including Aui ab tl do cancel identically such 
that the momentum-energy tensor is given in Eq. (2.7). 

Appendix B. Equations of motion 

We start from the Dirac equation 

-r (V + ^-A = mi> + 26 (T 5 T a ^)J a . (4.8) 
a \ 2a J 

The conjugate can be easily derived, but we do not present it here. Let us remember the 
definition of the pseudo scalar P = itpT^tp, scalar S = tpip, vector current V a = ipTaip, axial 
current J a = ■0r 5 r a ^, Lorentz tensor = §Y>[r o , Tf,]^, the factors of i are chosen such 
that they are all real. 

The first step is to establish the equation of motion for the scalar. In order to do so we 
multiply (4.8) by ipFo on the left and obtain 

- (W + = WW + 26 (^r r 5 r a v>) J a • (4.9) 

Since T§ = iToTiT2T3 we have that TqT^ = ^ril^r^ and Tor5r a = zril^rsEj = 
— ? ?aor5 — z/2e a f ? TT- J . Thus Eq. (4.9) becomes 

a + lr^) = mV ° + 26 ( iJ ° p ~ \ eai i Lijja ) ' (4 - 10) 

while its conjugate gives 

- l - (W + =mV + 26(- iJ°P - ^e ai] L^J a ) . (4.11) 

By taking the imaginary part of (4.10) we get the conservation law 

7^{a 3 S)' = 26J P. (4.12) 
By taking the real part of (4.10) we get 

± - V^V) = mV - 6 e aij L^J a . (4.13) 



- 12 - 



Now we look at the pseudo scalar, by contracting (4.8) with tpT^To we obtain: 
i ( 7 i 3a' 



Ut s i// + f^i^J = m^r 5 r v> + 29 (vT 5 r r 5 r a vO r . (4.u) 

Using r| = 1 and 'tpToFatp = ^rjoa^ ~~ *-^0a, the RHS reads m Jo — 29 (JoS — iLo a J a ^ . Taking 
the real part we obtain 

(« 3 ^)' = mJ -29J S. (4.15) 

We now consider the components of the axial vector by contracting (4.8) with ipT^ 
if- 3a' - \ 

- Ur 5 iW + — $r B iw J = m^r 5 t/> + 2d (^r 5 r 5 r a ^) r . (4.16) 

The RHS is —imP + 29V a J a , thus the imaginary part gives 

^4 (a 3 Jo)' = -mP. (4.17) 
We now look at the spatial components of the axial vector by contracting (4.8) with t^T^TiFo: 
% - UTsTiil/ + I^VTsr^) = m^TiToi; + 29 {^T 5 Y iT Q T 5 Y J a . (4.18) 

Since r 5 rjr = -i/2e ijk T j T k , and rjr r a = rj 0a Ti - rjiaJ^o + ie iab T 5 F b , the RHS reads 
-l/2me ijk V k + 29 (ViJ - JjVb + ie iab J a J b ). Thus taking the imaginary part we get: 

-(a 3 J;)' = 29e iab J a J b = 0. (4.19) 
a 

We now consider the time component vector current by contracting (4.8) with ip: 

% - (VoV' + f^ r oV>) = rrvH> + 29 (V5r 5 r a V) J a . (4.20) 
the RHS reads mS + 29 J a J a which is real, thus we obtain 

^ W = 0, (4.21) 

while 

- (ipTorp 1 - $'Toil>) =mS + 29 J a J a . (4.22) 
a 

The equation for the space component of vector current by contracting (4.8) with v/jTiTo: 

a + 1^^) = m ^ TlT ^ + 29 {^oTsTai,) J a . (4.23) 

Now rjr r 5 r a = -rjr r a r 5 = -rj 0a TiT 5 + i] ia r T 5 - ie aij T j , thus the RHS is -imL i0 + 
29iei a jJ a V :) . Thus we have the evolution equation 

^ (a 3 Vi)' = mL 0l + 29e iaj J a Vi (4.24) 
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while 

" (tpr^ 1 - ftr^) = o. (4.25) 



i 

a 



The equation for the space-time component of Lorentz tensor is obtained by contracting 
(4.8) with $Ti: 

if- V - \ 

- UViTo^' + —^TiToi)] = mijTii; + 29 {i>T i T 5 T a ^)J a . (4.26) 

Since Tir 5 r a = -i lia T 5 - ie mb T F b + i/2r] a0 e ijk T^T k the RHS is mVi + 29{iJiP - e iab L ob J a + 
l/2eijkL^ k Jq). Taking the real part of the equation we obtain 

^4 (« 3 W = mVi - 29e mb J a L ob + 9J e ijk V k . (4.27) 

The equation for the space components of the Lorentz tensor is obtained by contracting 
(4.8) with ^r^To: 

~ UviTjij)' + ^LjTiTji/A = mijT.TjTo^ + 26 (^^5^) J a . (4.28) 
Using that TiTjVQ = tj^Tq — ieij k T^r k and that 

rir^rorsra = -(rj ai ej bc - ri aj e ibc )T b r c - ir) a0 e ijk r°r k + ie ija - %a%T 5 - ^r] ij e ak iT k T l , 
we obtain, antisymmetrizing in i and j, 

^4 (a 3 ^i)' = 2e{^J t e jab L ah - \jje iab L ab - J e ijk L ok ) . (4.29) 



^j(a 3 L k Y = 29 (V' , ./,/. ; + J L ok ) . (4.30) 



Contracting this with e ijfc /2 gives 

1 

2^ 
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